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ABSTRACT 

We present numerical time-dependent calculations for fall-back disks relevant for GRBs in which the disk of 
material surrounding the black hole (BH) powering the GRB jet modulates the mass flow, and hence the strength 
of the jet. Given the initial existence of a small mass < 10 _4 M 0 near the progenitor with a circularization 
radius ~ 10 iO — 10 11 cm, an unavoidable consequence will be the formation of an “external disk” whose outer 
edge continually moves to larger radii due to angular momentum transport and lack of a confining torque. For 
long GRBs, if the mass distribution in the initial fail-back disk traces the progenitor envelope, then a radius 
~ 10 11 cm gives a time scale ~ 10 4 s for the X-ray plateau. For late times t > 10 7 s a steepening due to a 
cooling front in the disk may have observational support in GRB 060729. For short GRBs, one expects most 
of the mass initially to lie at small radii < 10 s cm; however the presence of even a trace amount ~ 1O _9 M 0 
of high angular material can give a brief plateau in the light curve. 


Subject headings: gamma rays: bursts 


1. INTRODUCTION 

The tentative pr e-Swift hints that breaks occur in the long- 
term afterglow light curves in different wavebands at ~ 1 — 10 
d (Frail et al. 2001) after the GRB have not been borne out 
by numerous detailed observations of Swift GRB afterglows in 
recent years (Ghisellini et al. 2007, Oates et al. 2007, Racusin 
et al. 2008, Liang et al. 2008). With Swift (Gehrels et al. 
2004), the long-term ( t > 10 6 s) X-ray behavior has turned 
out to be surprisingly complex, exhibiting alternating steep 
and shallow slopes when plotted as log F x versus log(£ — TO) 
(Zhang et al. 2006; Nousek et al. 2006). Recent workers have 
begun to explore the possibility that the long-term decay of 
the X-ray flux is not due to the deceleration of baryonic ejecta, 
but rather a secular decrease in the rate of accretion powering 
the central engine, and therefore indirectly the jet (Kumar et 
al. 2008a, 2008b, Metzger et al. 2008, Cannizzo & Gehrels 
2009, hereafter CG09). For long GRBs (lGRBs), the early, 
steep rate of decay may be giving us information about the 
radial density distribution within the progenitor core (Kumar 
et al. 2008b), whereas the later decay may be governed by 
the outward expansion of the transient disk formed from the 
remnants of the progenitor CG09). 

As regards short GRBs (sGRBs), even if a small amount of 
material (~ 10~ 5 - 1O" 4 M 0 ) is expelled during the NS-NS 
merger and later accreted in a disk, that would be sufficient 
to power a bright afterglow, which may also be strongly in- 
fluenced by the effects of r-process nucleosynthetic heating in 
the neutron rich material that becomes the disk (Metzger et al. 
2010 ). 

Zhang et al. (2006) present a schematic for the decaying 
GRB light curve as seen by the XRT on Swift. The decay 
is traditionally shown in logF - logt There are four basic 
power-law decay (F ex t~ a ) regimes: (i) a steep decline fol- 
lowing the prompt emission with ai cz 3 out to 10 2 — 10 3 s, 
(ii) a plateau with a n — 0.5 out to 10 3 — 10 4 s, (iii) a steep- 
ening with am ~ 1.2 out to 10 4 - 10 5 s, and (iv) a further 
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steepening at late times (not always seen) with a\y — 2. 

CG09 present a general analytical formalism to explain the 
different power law decays using a fall-back disk, where the 
variations in a could potentially be accounted for by different 
physics operating within the disk. The results of CG09 were 
purely analytical; in this work we present time dependent nu- 
merical calculations in order to examine in more detail the 
potential of the model, and we apply the results to XRT data 
for one 1GRB and one sGRB, taking the best studied of each 
class. In Section 2 we review the Dainotti relation, an empiri- 
cal relation involving the duration and luminosity of segment 
II, in Section 3 we present our detailed numerical model, in 
Section 4 we compare the model with observations for the 
1GRB 060729, the GRB with the longest observational time 
series in X-rays, in Section 5 we compare theory and data 
for the sGRB 05 1221 A, the sGRB with the longest and most 
detailed XRT light curve, in Section 6 we revisit the Dainotti 
relation, in the context of our numerical results, and in Section 
7 we discuss and summarize our results. 

2. THE DAINOTTI RELATION 

Dainotti et al. (2008, 2010) found an empirical relation be- 
tween the duration of the X-ray plateau in the source frame 
= tn(l + z)"" 1 , and the X-ray luminosity at the end 
of the plateau, also corrected into the source frame, Ljj. 
Expressing their relation in the form log 10 Tjj(erg s" 1 ) = 
a -j- 61og 10 tjj(s), Dainotti et al. (2010) find a = 51.1 ± 1 
and b = — 1.1 ± 0.3. Using the same sample of 62 lGRBs 
with small errors taken from Dainotti et al. (2010) we find 4 
a = 50.5 ± 0.4 and b = —0.92 dt 0.1, close to their results. 

If the inverse relation between L and is physical, and 
if long-term accretion is the correct explanation for the long- 
term X-ray light curve of GRBs, that would have implications 
for the amount of mass in the initial fail-back disk. To first 
order, the fact that L^t jj is constant would imply the accreted 
mass reservoir is constant. Willingale et al. (2007) present a 

4 We employ a different method than Dainotti et al. (2010), a Monte Carlo 
technique in which 10 6 data sets are created with ler errors randomly either 
added or subtracted to each data point in both x and y directions. Thus for N 
data points one could in principle have N 4 distinct data sets. For each data 
set, a and 6 values are calculated, and the final values for a, <5a, b, and Sb are 
taken from the averages and standard deviations of the 10 6 values. 



2 


CANNIZZO, TROJA, & GEHRELS 



Fig. 1 . — Inferred tStal accretSn mass fof the platea% + later cficay phases 
of GRBs, using data from Dainotti et al. (2D 10) for their 62 well-constrained 
lGRBs. The hatched area (shown in red) indicates a putative limiting XRT 
detection flux level fu ~ 10" 12 erg cm” 1 s" 1 (adopting a plateau duration 
£jj = 10 4 s) for being able to study a plateau to sufficient accuracy that it 
would have been included in the Dainotti et al. sample of 62 GRBs with 
good statistics (i.e., u = [(£log£jj) 2 + (^logfjj) 2 ] 1 / 2 < 4) and known 
redshifts. We adopt a beaming factor / == 1/300 and an accretion efficiency 
for powering the X-ray flux e acc — 0.03 to convert from X-ray fluence (i.e., 
total X-ray energy, after taking into account 47rri 2 ) to accreted mass. 

simple formalism for integrating the X-ray flux on the plateau 
and subsequent am ~ 1.3 decay to obtain a total energy. 
We can convert this into a mass by making a few plausible 
assumptions about the energetic efficiencies (Krolik, Haw- 
ley, & Hirose 2007), e.g., an X-ray afterglow beaming factor 
/ ~ 10~ 3 — 10“ 2 , and an efficiency e aC c — 0.01 — 0.1 with 
which accretion onto the inner engine powers the observed X- 
rays, presumed to be created within the beamed jet. We use 
eqns. [3] and [4]) from Willingale et al. (2007) to obtain a 
total energy 5 Ex and hence 5M derived from the X-ray flu- 
ence. The results are shown in Figure 1. The red hatched 
area indicates The hatched area is bounded from above by an 
estimate of the effective limiting XRT flux 5 for being able to 
observe and characterize plateaus sufficiently well that they 
would satisfy u < 4. There appears to be a selection effect 
giving the Dainotti relation, namely the flux detection limit 
for XRT prevents one from observing faint plateaus at high z. 
The Dainotti relation, effectively equivalent to ~ con- 
stant, results from the sampling of GRBs primarily beyond 
z ~ 1.5, for which faint plateaus are observationally biased 
against. Thus we are seeing basically a narrow strip corre- 
sponding to the upper end of a much broader distribution, 

5 Our adopted XRT plateau detection level 10" 12 erg cm" 1 s" 1 signif- 

icantly exceeds the nominal single pointing detection limit ~ 10" 14 erg 

cm" 1 s" 1 because one must take into account the effect of Swift orbital 

gaps, and also the need to distinguish the different phases of the canonical 

light curve (i.e., plateau versus steep decay). In addition, the presence of 

flares may further complicate the picture, but they usually are confined to 

t < 10 3 s and generally do not represent a major portion of the total energy 

budget. 


which makes it appear that ~ constant. For z < 1.5 

one sees curvature in the lower envelope of SEx values due 
to the strong dependence of the detection limit on z. Never- 
theless, the fact that there appears to be a well-defined upper 
limit SM « 10“ 4 — 10“ 3 M® is interesting: if one starts with 
a ~ 100 Mq progenitor and if accretion governs the long-term 
X-ray light curve, this indicates that no more than a fraction 
~ 10~ 6 — 10 ~ 5 of the progenitor mass survives the hypemova 
to form a fail-back disk. 


3. ACCRETION DISK PHYSICS 

By writing the equations for mass continuity and angular 
momentum transport in cylindrical coordinates, assuming Ke- 
plerian rotation h 2 K — GMbht” 3 and integrating over the 
vertical thickness of the accretion disk, one arrives at an equa- 
tion for the evolution of the surface density X = 2 ph, where p 
is the density and h the disk semithickness (actually pressure 
scale height), 


(^Er 1 / 2 ) 


The kinematic viscosity coefficient 

2aP 

K p' (2) 

where P is the pressure and a is the Shakura-Sunyaev 
parametrization of the angular momentum transport and heat- 
ing (Shakura & Sunyaev 1973). Equation (1) is discretized 
following the method of Bath Sc Pringle (1981), in which grid 
points are distributed as r 1 / 2 . In addition, one has a thermal 
energy equation governing the temperature evolution, 

dT 2 (A-B + C + D) JIT Id, . . dT ^ 

-57 = — (ru r )-u r — , (3) 

at c p H pc p r or or 

where the viscous heating A = (9/8) ^Q 2 S, the radiative 
cooling B — <jT * , and C and D represent radial heat fluxes 
due to turbulent and radiative transport (Cannizzo et al. 201 0). 
For the calculations we present in this work, the thermal equa- 
tion only becomes of potential importance at very late times 
> 1 yr. 

CG09 present analytical models for the afterglow light 
curves powered by fall-back disks. The main strength of 
model is the universal decay law dlogL &cc /dlogt ~ —1.3 
characterizing fall-back disks without the external, confining 
tidal torque of a companion star (Cannizzo et al. 1990). This 
law seems consistent with the late time decay of GRBs seen 
in X-rays. The early time decay, in particular the steep decay 
for t < 10 3 s, following the prompt emission, is difficult to 
explain within the fail-back disk scenario and may be giving 
us information about the radial density profile within the pro- 
genitor (Kumar et al. 2008b). In this work we present results 
of time dependent numerical calculations using a general ac- 
cretion disk code (Cannizzo et al. 2010), which uses the input 
physics detailed in CG09. The boundary conditions are as fol- 
lows: (i) at the inner disk edge, taken to be 10 7 cm for sGRBs 
and 10 8 * * * cm for lGRBs, matter is instantly removed as it ar- 
rives, while (ii) the outer grid point is placed at such a large 
radius, typically 10 13 cm, that during the course of the run 
the outer edge of the spreading fail-back disk never reaches it. 
The large dynamic range in disk radii r outev /ri nneT necessi- 
tates ~ 1500 — 3000 grid points. In addition, no fresh material 
is added during a run. Thus the evolution is set entirely by gra- 
dients within the disk, unlike the standard Shakura-Sunyaev 
disk fed at a constant rate in the outer edge which approaches 
a steady-state M(r) = constant with time. 
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Fig. 2. — Model lig^cmj*ve showing the accretion-defrf$P1&-ray flux 
from the fall-back disk in a 1GRB, taking Mbh = ISM© and Mdisk = 
10 " 4 M©. Also shown is the Swift XRT light curve (Evans et al. 2007, 
2009) of GRB 060729 (in blue). Given the value 4 tt d\ = 1.15 x 1G 57 
cm 2 for GRB 060729 ( z — 0.54), the efficiencies assumed in the two scaled 
model light curves (in red) to convert from rate of accretion onto the central 
engine to XRT fluxes are ft = 8.8 x IQ" 3 (upper) and 1.4 x 10" 2 
(lower). The five flux values at late times 10 7 s < t < 10 s s (in green) come 
from Chandra observations (Grupe et al. 2010). The last data point lies at 
t — TO = 642 d. The flat decay portion in the model light curve corresponds 
to the time during which the initial £(r) is being redistributed into an accre- 
tion disk with a self-consistent radial profile. The bottom panels show £(r) 
(solid lines, top to bottom) at t = 5, 25, 100, 300, and 600 d. The dashed line 
indicates the initial £(r) distribution, and the parallel dotted lines indicate 
the critical values associated with the dwarf nova limit cycle which occurs at 
the transition between neutral and ionized gas. 

4. LGRBS 

LGRBs are thought to be due to the explosion of a short- 
lived > 30 M© progenitor (MacFadyen & Woosley 1999). 
The hosts for lGRBs tend to be subluminous, irregular galax- 
ies rich in star formation (Fruchter, et al. 2006). If the long- 
term X-ray light curves of GRBs are indicative of feeding 
from a fall-back accretion disk, then from Figure 1 we see 
that no more than ~ 10“ 4 - 1CT 3 M© of material in the fall- 
back disk is needed, for nominal assumptions about the effi- 
ciencies. In other words, for a ~ 100M© progenitor, a mass 
fraction only ~ 10~ 6 - 10 -5 (excluding the ~ 10M© which 
ends up within t < 10 - 10 2 s in the BH) is required to sur- 
vive and persist in the vicinity of the progenitor to power the 
long-term light curve. 

Figure 2 shows a light curve for Mbh = 15M© and 
Mdisk = 10~ 4 M©. The initial E(r) profile is taken to be 
ex r“ 2 / 3 out to 7 x 10 10 cm, and zero for larger radii. The 
lower panels reveal the £(r, t) evolution. The steep decay as- 




Fig. 3. — Model lighr(xsm) showing the accretion-derivedT^ppJwer from 
the fail-back disk in a sGRB, taking Mbh — 3M© and Mdisk — 10" 5 M©. 
Also shown is the Swift XRT light curve (Evans et al. 2007, 2009) of GRB 
05 1221 A (in blue). Given the value 4tt d\ = 1.18 x 10 57 cm 2 for GRB 
05 1221 A (z = 0.5465), the efficiencies assumed in the two scaled model 
light curves (in red) to convert from rate of accretion onto the central engine 
to XRT fluxes are ftacc = 2.4 x 10" 2 (upper) and 4.3 x 10" 2 (lower). One 
sees fairly rapid decay initially due to the presence of a significant amount of 
accreting material at small radii. The brief shoulder in the light curve is due 
to the delayed accretion of the spiral arm of ejected NS matter at large radii, 
motivated by detailed SPH calculations (Rosswog 2007). The bottom panels 
show the surface density evolution, at equally spaced intervals St = 2d, The 
dashed line indicates the initial £(r) distribution, and the parallel dotted lines 
indicate the critical values associated with the dwarf nova limit cycle. 


sociated with segment I is not accounted for in our model, and 
may well be due to a separate physical process, such as the ac- 
cretion of the progenitor core (Kumar et al. 2008b; Lindner et 
al. 2010, see their Fig. 2). The time scale tn for the plateau 
associated with segment II comes out naturally in the models, 
given the ~ 10 11 cm radius for the progenitor — assuming that 
the initial fall -back mass distribution roughly traces the enve- 
lope of the progenitor. The observed X-ray decay for GRB 
060729, the GRB which has the been observed for the longest 
time in X-rays (Grupe et al, 2007, 2010), does not match pre- 
cisely that of the fall-back disk. There may be various system- 
atic effects that could account for the difference. The most ob- 
vious effect, a variable cosmological K-correction concomi- 
tant with the fading in X-ray flux over a dynamic range of 
~ 6 decades, cannot play a strong role, as the spectral index 
does not vary significantly over the long term from its mean 
value 0 ~ 2. There may be some time -variable physics asso- 
ciated with the X-ray emission from within the jet that could 
affect either e acc or / versus time. Our calculations only rep- 
resent the accretion disk which powers the jet and we take 
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6 a cc = / = 0 for simplicity. Modeling the jet emission is 
beyond the scope of this work. 

At late times 10 6 7 s < t < 10 8 s there is a deviation from 
the canonical power law decay because S(r ou ter) drops to the 
point at which the dwarf nova limit cycle becomes active (in- 
dicated by the parallel dashed lines), and a cooling front is 
launched from the outer edge that instigates a transition from 
ionized to neutral gas. This quenches the source of accretion 
onto the central engine. The very late time Chandra observa- 
tions of GRB 060729 (Grupe et al. 2010) appear to coincide 
with the timing of this drop-off, which may therefore repre- 
sent a cooling transition front in the outer disk which dimin- 
ishes the rate of accretion onto the central engine, rather than 
a jet break. 

5. SGRBS 

SGRBs are thought to be caused by the merger of two neu- 
tron stars (Eichler et al. 1989; Paczynski 1991; Narayan, Pi- 
ran, & Kumar 2001; Rosswog & Ramirez-Ruiz 2002). The 
hosts for sGRBs tend to be ~ 5 times more spatially extended 
than those for IGRBs (matching the ~ 5 times wider projected 
spatial distribution of sGRBs vs. IGRBs within their hosts), 
are also more indicative of an older population with less ac- 
tive star formation (Fong, Berger, Sc Fox 2010). Due to their 
faintness relative to IGRBs, sGRBs are much less well-studied 
in X-rays. Indeed, there is only one certifiable example of a 
sGRB which has a long and well-sampled X-ray light curve, 
to the extent that a statement can be made as to the presence 
of a plateau - GRB 05 1221 A. 

The importance of the outward spreading of the disk formed 
from the NS-NS merger was clearly demonstrated by Met- 
zger et al. (2008, see their Fig. 1). The initial radial density 
profile, and hence X(r) profile after after angular momentum 
conservation has vertically compacted the merger remnants 
into their accretion plane, is of course much more compact 
radially than that expected for IGRBs. Theoretical guidance 
on a potential E(r) profile to begin the calculations, particu- 
larly at large radii, is scarce, yet there are indications of small 
amounts of matter ~ 10~ 6 — 10~ 4 M© ejected at early times 
that may have high angular momentum, and therefore circu- 
larize at large radii (Rosswog 2007). 

Figure 3 shows the evolution of a fall-back disk of poten- 
tial relevance for the aftermath of a NS-NS merger. We take 
M B h = 3M© and Mdisk = 10“ 5 M©. Unlike the much 
longer evolution shown in Figure 2 for IGRBs, in this case 
the outer edge of the disk is still freely expanding to larger 
radii by the end of the run (indicated by the narrow spike in 
£(r) between 10 11 cm and 3 x 10 11 cm shown in the second 
lower panel). The brief plateau at ~ 10 4 s results from the 
ad hoc introduction of the small amount ~ 10 ” 9 M© of high 
angular momentum material at the large circularization radius 
10 10 cm. The high efficiencies associated with BH accretion 
show the potential for a small amount of material < lAf© to 
have a dramatic effect on the long-term light curves as regards 
brief plateaus or inflection points. 

6 . POTENTIAL APPLICATION TO A DAINOTT1-LIKE RELATION 

In Section 2 we show that the Dainotti relation L\ \ 

as originally envisioned may be due to the observational bias 

against detecting and characterizing faint plateaus. Neverthe- 
less, the accreted mass estimates SM inferred using the LJj 

and t{ j values from Dainotti et al. (2010) are interesting in 

the context of this work, and one might legitimately ask what 



Fig. 4. — Model light curves for 1GRB parameters, keeping the initial fall- 
back disk mass constant at 1 O" 4 M 0 but varying the initial radius and overall 
normalization. For the central light curve ro = 7 x 10 10 cm, as in Figure 
L For each successive light curve ro is increased a factor of two going to the 
right, and decreased a factor of two going to the left. To keep SM constant, 
the overall normalization on E(0) is varied by an additional factor of 2.5 for 
each successive run. 

theoretical prediction for the Ljj(£ Jj) relation the fall-back ac- 
cretion hypothesis would make, given a hypothetical physical 
constraint in which SM were held constant and the initial ra- 
dius of the fail-back disk were allowed to vary. Figure 4 shows 
the results of seven runs in which we fix SM = 10 “ 4 M© and 
vary the initial outer radius of the fall back disk. One can see 
a clear inverse relation between the duration of the plateau j 
and the luminosity at the end of the plateau L\ v The low-;? 
behavior of SM seen in Figure 1, however, indicates that a 
spread in SM may be more realistic, in which case one would 
not expect to be able to use the theoretical prediction of the 
Tjj(£jj) relation as a useful discriminant for the theory. 

7. DISCUSSION AND CONCLUSION 

We have presented time dependent calculations of the fall- 
back disk scenario to account for the long-term X-ray light 
curves for GRBs. For IGRBs, an initial radial scale of a stel- 
lar radius ~ 10 11 gives a natural viscous evolution time of 
~ 10 4 s to redistribute the matter into a quasi-steady disk, 
roughly consistent with the observed plateau duration. The 
rate of decay for GRB 060729 is close but does not match 
the models precisely, which may hint at time-variable emis- 
sion processes that would affect our adopted efficiencies e acc 
and /. It is interesting to note that our effective values of the 
efficiencies, accretion plus beaming, required to match the ob- 
served X-ray flux levels are comparable between IGRBs and 
sGRBs: /e™* - KT 2 for GRB 060729 and fe^ c ~ 3x 10” 2 
for GRB 051221 A. If the accretion efficiency e acc is about the 
same for IGRBs and sGRBs, this may indicate that sGRBs are 
less beamed by a factor of ~ 3 compared to IGRBs, roughly 
in line with previous results (Watson et al. 2006, Grupe et al. 
2006). This similarity between IGRBs and sGRBs afterglows 
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is not entirely unexpected, given the general similarities in 
their afterglow properties (Gehrels et al. 2008; Nysewander, 
Fruchter, & Pe’er 2009). 

At very late times Chandra observations indicate a steepen- 
ing in the rate of decay (Grupe et al. 2010), which appears to 
be consistent with the onset of a cooling front in the disk. This 
would be an alternative to the standard jet-break interpretation 
discussed by Grupe et al. (2010). For sGRBs, the picture is 
less clear, given that we only have a single well-studied ex- 
ample. We have shown that the presence of even a very small 
amount of high angular momentum gas ~ 10“ 9 M© can give 
a slight inflection to the X-ray decay, as was observed in GRB 
05 1221 A. If only ~ 10 _4 M© of gas survives either the hy- 
pemova (lGRBs) or NS-NS merger (sGRBs), then the accre- 
tion resulting from the ensuing fail-back disk should power a 
long-term jet. 

Lastly, we have shown that the Dainotti relation Ljj oc ^n 1 
may be due to an observational bias against detecting and 
characterizing faint plateaus: the relation is governed by 
GRBs at 2 > 1,5 for which we only detect the upper enve- 
lope of a broad distribution. Nevertheless, the existence of 
an apparent upper limit to the total X-ray energies inferred 
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if one assumes accretion onto the central engine as the long 
term powerhouse for the X-ray flux, is extremely interesting. 
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mum mass fraction ~ 1Q~ 6 — 10” 5 of the progenitor survives 
in the vicinity of the progenitor to be accreted as a fail-back 
disk (excluding the ~ 10M© that ends up in the BH during 
the prompt emission and subsequent segment I consisting of 
the steep-decay). This has important ramifications for the en- 
ergetics associated with the hypernova explosion and subse- 
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